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ABSTRACT 

Recent observations with Galaxy Evolution Explorer (GALEX) show strong unexpected UV excess in the 
spectrum of brightest cluster galaxies (BCGs). It is believed that the excess UV signal is produced by old and 
evolved core-He burning stars, and the UV flux strength could be greatly enhanced if the progenitor stars have 
high value of He abundance. In this work, we propose that sedimentation process can greatly enhance the He 
abundance in BCGs. Our model predicts that the UV flux strength is stronger in more massive, low-redshift, 
and dynamically relaxed BCGs. These predictions are testable with the current generation of GALEX+SDSS 
observations. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: clusters: general — ultraviolet: galaxies 
— diffusion 



1. INTRODUCTION 

Brightest Cluster Galaxies (BCGs) are among the most 
massive galaxies in the Universe, forming and evolv- 
ing at centers of galaxy clusters through galaxy mergers 
and accretion. The baryonic component of present-day 
BCGs is composed of mainly old stellar p opulations (e.g. 
iDe Lucia & Blaizotll2007l: ICollins et al.ll2009l) with a fraction 
of younger stellar populations formed out o f cooling gas in 
cluster cores (e.g.. iPeterson & Fabian! 12006). Properties of 
BCGs are therefore dependent on the mass accretion histories 
and the conversion of hot gas in clusters into stars. Detailed 
studies of the BCGs properties have the potential to shed new 
insights into the physics of the most massive galaxies in the 
Universe. 

The UV upturn is one of the mysterious phenomena ob- 
served in massive elliptical galaxies. Since the first dis- 
covery ( ICode & Welchlll979l) . it has been known for several 
decades that some giant elliptical galaxies exhibit strong un- 
exp ected bump in the Ultraviolet (UV) part of the spectrum 
(see IO'Connelll[l999l for a review). Recent GALEX obser- 
vations in dicate that the U V flux excess is most significant 
in BCGs dRee et al.ll2007b . The observed UV excess, also 
know n as the UV upturn, is spatially smooth and extended 
(e.g.. lLee et afll2005at IRee et al.ll2007l) . The spatial distribu- 
tion and the spectral shape of the UV upt urn phenomena can - 
not be explained by recent star formation dBrown et al.ll 1 997b . 

One of the leading models of the UV upturn suggests 
that the phenomena can be explained most naturally by old 
evolved stars, known as extreme horizontal-branch (EHB) 
stars. EHB stars are hot, core helium-burning stars with ex- 
tremely thin hydrogen envelopes (M env < 0.05M o ), which lost 
their envelopes via massive winds o n the red giant branc h in 
the single star evolution model (e.g.. lDorman et af]|1993l) 3 . 

High value of He abun dance can make the formation of ho t 
EHB stars more effective dDorman et al.ll995UYi et al.ll997l) . 
This is because that HB stars with a higher helium abundance 
evolve faster on main sequence and red giant branch, can burn 
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more hydrogen during the core helium-burning phase, and 
therefore can become UV-bright more easily. In particular, 
a very large He abu ndance (Y > OA) c an greatly boost the 
UV upturn strength dDorman et alJI 1 99~5h . 

In this work, we point out that helium sedimentation pro- 
cess can greatly increase the helium abundance in BCGs, 
making them conducive to the formation of UV-bright EHB 
stars. It has long been suggested that helium sedimenta- 
tion occurs in the intracluster medium (ICM) of galaxy clus- 
ters (lAbramopoulos et al.l 11981b iGilfanov & Svunvaevll 19841: 
| Oin & Wull2000UChuzhov & Nusseril2003HChuzhov & Loebl 
l2004t lEttori & Fabianll2006tlPeng & Nagaill2009l) . Under the 
influence of gravity, the heavier helium nuclei in the H-He 
dominated ICM accumulate at the center of massive galaxy 
clusters. If the UV flux is produced mainly by EHB stars, the 
He sedimentation process can boost the UV flux produced by 
these stars. Our sedimentation scenario predicts that the UV 
upturn phenomena should be most pronounced in high-mass 
(or Tx), low-redshift, and dynamically relaxed systems. We 
show that these predictions are testable with the current gen- 
eration of GALEX+SDSS observations. 

2. HE SEDIMENTATION EST GALAXY CLUSTERS AND ELLIPTICAL 
GALAXIES 

2.1. Diffusion Equations 

We calculate the evolution of the He sedimentation in 
galaxy clusters and elliptical galaxies by solving the diffusion 
equations for a fully ionized hydrogen and helium plasma. 
Detailed descriptions of these calculations are described in 
dPeng & Nag ai 2009). Here we provide a brief overview of 
the calculations. 

Each species s obeys an equation of continuity and momen- 
tum conservation, 

dn s | 1 d(r 2 n s u s ) _ Q 
dt r 2 dr 
dP s v-i 

— + n s A s tn p g - n s Z s eE = ^ K st {w, - w s ) . (2) 

Here the species s has mass A s m p , charge Z s e, density n s , 
partial pressure P s , and velocity u s , where m p is the proton 
mass. The effective resistance coefficients is defined as K st = 
f B l Kf t =a , where fs is a magnetic or turbulent suppression fac- 
tor and Kf. =0 oc j^l 1 is the resistance coefficient of the un- 
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magnetized H-He plasma ( IChapman & Cowlindl952l) . In the 
following calculations, we assumed that — 1 to obtain the 
maximum He sedimentation effect. Equation (f2]i describes 
forces acting on a specie s, and it is the balance of these forces 
that ultimately determines the rate of sedimentation. For a 
sinking He nucleus, the gravitational force (g) is counteracted 
by three types of forces provided by the induced electric field 
(E), the pressure gradient (dPJdr) of helium, and the drag 
force due to collisions with surrounding particles. Note that 
the sedimentation destroys hydrostatic equilibrium since re- 
distribution of particles introduces a temporal change in the 
total gas pressure. However, hydrostatic equilibrium can be 
restored quickly. This equilibrium restoring acquires a net in- 
flow with a mean velocity u = £ s n s A s u s / Yj S n sA s , 

du 1 <9P gas 

-77 = o g, (3) 

dt Pga S or 

where P gas = Yj S n sksT is the total gas pressure for ideal gas, 
and p eas = 2.s n s A s m u is the gas density. 

The diffusion velocity between species s and the fluid ele- 
ment is w s — u s - u, which satisfies mass and charge conser- 
vation, 

Y j A s n s w s = 0, (4) 

S 

^Z s n s w s =0. (5) 

Note that the summations include both ions and electrons. To 
satisfy these conservation laws, for each sinking helium nu- 
clei, there are roughly four protons and two electrons that float 
up. Equations (Q]i-([5]i describe the process of particle diffusion 
in a multi-species plasma. 

2.2. Cluster and Galaxy Models 

We set up cluster and galaxy models and initial conditions 
as follows. Initially, we assume that the gas consists of a pri- 
mordial H and He plasma uniformly throughout clusters (or 
galaxies) (X = 0.75 and Y = 0.25). We ignore the contribu- 
tion of elements heavier than He. We set up the initial gas dis- 
tribution by assuming hydrostatic equilibrium of the plasma 
medium in the potential well of clusters dominated by dark 
matter, following the Navarro-Frenk- White (NFW) densit y 
profile for the total mass distribution (Na varro et alj|1997h . 
We adopt t he concentration para meter C500 = rsoo/r. s = 4 
for cluste rs (IVikhlinin et aT]l2006l) ("V06) and 10 for elliptical 
galaxies (IWechsler et al.l 120021) . where r s is the scale radius 
of the NFW density profile and rsoo is a radius of a spherical 
region within which the mean enclosed mass density is 500 
times the critical density of the universe. 

For the temperature profile, we consider the observed tem- 
perature profile obtained with deep Chandra observations of 
nearby relaxed clusters (V06). The observed temperature 
peaks around r/rsoo - 0.2, and decreases at both inner and 
outer radii (see Figure 1 in iPeng&Nagail 12009). The ob- 
served temperature drop in the inner radii leads to a signifi- 
cant suppression of He sedimen tation, compared to the r esults 
based on the isothermal model (IChuzhoy & Loebll2004l) . For 
elliptical galaxies, due to their lower X-ray luminosity, the 
temperature observations are more difficult and usually lim- 
ited to core regions. We assume that they have the same uni- 
versal temperature profiles of clusters. Most of the observed 
brigh t elliptical galaxies do show V06-like tempe rature pro- 
files dDiehl & Statlerll2008t IFukazawa et al.l 120061) . The gas 



density is derived from solving the hydrostatic equilibrium 
equation. The central density is normalized by requiring the 
enclosed gas mass fraction to be / ga s(< rsoo) = 0.15 for clus- 
ters (V06) and 0.015 for elliptical galaxies. 4 Throughout this 
work, we use cosmologica l parameters: Qm = -3, Q.\ = 0.7, 
Q b = 0.0462, and h = 0.7 dKomatsu et al.ll2008l) . 

2.3. Sedimentation Velocity and Scaling Relations 

To develop physical insights into the process of He sedi- 
mentation, it is useful to consider a drift velocity of a trace 
He particle in a background of hydrogen, i.e, « p » « H e and 
w p = 0. In this limiting case, equations (fJJ for the two species 
are decoupled. The right-hand side of the equations for H 
vanishes, thereby fixing an electric field, eE = 0.5 m p g. Sub- 
stituting E into the equation of motion for He, we obtain the 
sedimentation velocity of He nuclei as wee = 3 m p gnn e /K p u e , 
which gives 

WHe , 80 km s"> (_*_f * f * \ ( * V 
UOkeV/ \10- 7 - 5 cms- 2 /\10- 3 cm- 3 / 

(6) 

where the induced electric field counteracts gravity and sup- 
presses the sedimentation speed by 25%. At a fixed density, 
the sedimentation velocity is generally larger for higher tem- 
perature (T) and gravity (g). Note that the pressure gradient 
(dP/dr) term in equation (f2l) further suppresses the sedimen- 
tation process. 

Substituting Eq. (O into Eq. ([TJ and integrat- 
ing over time, we obtain Y - Fo^ 4 T > where 
To is the primordial He mass fraction and r 
10Gvr(^ s /0.1)(r/10keVT 3/ 2 (l + 1.5 din T/d ln r)" 1 
is the sedimentation timescale dChuzhov & Lo eb 2004). Note 
that the gas temperature gradient term, d In T/d In r, takes 
into account of the non-isothermal temperature profile. For 
t « t, the He enhancement factor, AT = Y — To, has the 
following scaling AT <x At/r oc AtfZ^T 3 ' 2 . This indicates 
that the He sedimentation process is efficient for systems with 
higher T and lower / gas . 

2.4. He Abundance Enhancement in BCGs 

In this section, we solve the diffusion process in the H-He 
plasma numerically. Figure[T]shows the spatial distribution of 
the He enhancement factor, AT, for a Tx = 10 keV galaxy 
cluster and a Tx = 1 keV elliptical galaxy. The figure shows 
that the He enhancement could be significant (AT ~ 0.3 - 0.4) 
in the core of massive, relaxed clusters. For clusters, the He 
enhancement factors peak around r 0.05 rsoo and reaches 
0.4 (0.2) for cluster age of 1 1 (5) Gyrs. For elliptical galax- 
ies, the He enhancement is also non-negligible (AT ~ 0.1), 
because of their low gas mass fraction. Here, the He enhance- 
ment factor peaks around 0.1 rsoo with AT = 0.1 (0.05) after 
1 1 (5) Gyrs. Because of their higher concentration parameter, 
the elliptical galaxy has significant suppression of He sedi- 
mentation in the galaxy central core region. This is because 
that the higher pressure gradient causes more efficient out- 
ward He diffusion and hence suppresses the process of He 
sedimentation. 

Figure [2] shows the He enhancement at r = 0.1 rsoo as a 
function of X-ray temperature of galaxy clusters (top panel) 
and elliptical galaxies (bottom panel). Note that 0.1 rsoo 
roughly corresponds to the effective (half-light) radius of stel- 

4 The baryon components of elliptical galaxies are dominated by stars. 
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Fig. 1. — The spatial distribution of He enhancement factor in a Tx = 10 
keV cluster {top panel) and in a Tx = 1 keV elliptical galaxy {bottom panel). 
The set of curves in each panel, starting from the flat line, are for ages of 
0, 1, 3, 5, 7, 9 and 11 Gyr, respectively. For galaxy clusters, we assumed 
/eas(< '"500) = 0.15 and C500 = 4. For elliptical galaxies, we assumed / e as(< 
fsoo) = 0.015 and C500 = 10. 

lar surface brightness for BCGs. 5 Figure [2] demonstrates that 
the He sedimentation is strongly sensitive to the X-ray tem- 
perature. The He enhancement is also inversely proportional 
to the gas mass fraction, leaving the He sedimentation in el- 
liptical galaxies finite. 

Since the BCGs are galaxies located at the center of clus- 
ters, the He enhancement in BCGs is roughly given by the 
sum of the effects in clusters and elliptical galaxies: AY ~ 
0.5(0.25) after 1 1(5) Gyrs. Note that this is the optimal He en- 
hancement in BCGs since the He sedimentation process could 
be suppressed by magnetic field and turbulence. 

3. IMPLICATION FOR THE UV PHENOMENA 



In § 12.41 we showed that He sedimentation can greatly en- 
hance the He abundance in BCGs and that the effect is very 
sensitive to the X-ray temperature of host galaxy clusters. 
As we discussed in Section Q] high He abundance can make 
the environ ment conducive to th e formation of the UV-bright 
EHB stars dDorman et al.ll 19951) . Thus, the He sedimentation 
may provide an attractive mechanism to explain why the UV 
upturn phenomena is pronounced in BCGs. 

Connecting the idea of the He sedimentation to observa- 
tions of the UV upturn phenomena is a nontrivial step that 
requires detailed modeling of star formation and chemical en- 
richment histories as well as the full stellar population syn- 
thesis calculation that take into account the EHB formation 



Fig. 2. — The dependence of the He mass fraction at r = 0.1 rsoo on X- 
ray temperature of galaxy clusters (top panel) and elliptical galaxies {bottom 
panel). The curves starting from the flat line are for ages of 0, 1, 3, 5, 7, 9 
and 1 1 Gyr, respectively. 



and the He sedimentation. The full self-consistent modeling 
is thus beyond the scope of this work. However, it would be 
useful to check if the He sedimentation can have prominent 
effect on boosting the UV flux in BCGs, more specifically, 
the observed Far-UV (FUV)-to-V flux ratio, also known as 
the UV flux strength. The FUV flux is contributed mainly 
by stars in the HB/post-HB phases, while the V flux is con- 
tributed mainly by main-sequence stars. The FUV-to-V flux 
ratio is therefore given by L FUV /L V = n(t)efuv, where the 
total FUV flux is governed by the rate of stars entering the 
HB phases (or leaving off the main sequence) («(0) 6 an d the 
energy output in FUV band over the lifetime of HB/post-HB 
evolutionary track (ffw). 

High He a bundance can greatl y enhance the UV energy 
outputs epuv dDorman et al.ll 1 995b . An increase of He abun- 
dance yields a smaller main-sequence turnoff mass at a given 
age, and therefore a smaller envelope mass. It was demon- 
strated that sufficiently He-enhanced younger stellar popula- 
tion can still have smaller Main Sequence (MS) turnoff mass 
than an older populatio n with a normal He abundance (see 
e.g jDorman et alJl995L Table 4). A larger He abundance can 
also cause stars with an given envelope mass to burn more of 
their hydrogen envelope during their He core burning phase. 
Helium-rich stars can therefore produce UV flux for a larger 
range of initial envelope masses. For example, increasing Y 
from 0.27 to 0.47 roughly increases the total FUV energy out- 
put by a factor of 4 for a uniform distribution of envelope 



3 For BCGs of ste llar mass 10^ M a , .1 r 500 
radii - 30- 100 kpc (Graham et aU2006l) . 



50 kpc and the effective 



' The number of stars leaving main sequence per unit time per unit Ly. 



4 



Peng & Nagai 



mass ( iDorman et al.lll99~5l) . h(t), on the other hand, is fairly 
insensitive to the compositions. For the range of He abun- 
dance considered in this work, the enhancement factor of the 
UV flux strength due to high He abundance is < 4. 

Assuming that most of the stars are formed in a star burst 
at an early redshift (denoted as 'org'), and the rest is formed 
continuously with a constant star formation rate out of He- 
enriched gas (denoted as 'con'), we can express the total FUV- 
to-V flux ratio as 




+ /cc 



(1 + a/con)- 



(7) 
(8) 



where / con is the ratio of the continuous forming stellar mass 
to the total stellar mass, and a is the enhancement factor of 
UV flux strength due to high He abundance. In deriving Eq. 
(0, we assumed that the visual flux is not significantly af- 
fected by the He abundance. This is reasonable since the 
young He-enhanced stellar population has similar MS turnoff 
mass as the old first-generation stellar population. 

While it is widely believed that AGN feedback may sig- 
nificantly su£2ress_j^rj^o^ 

galaxies dPipino et alJl2009at ISchawinski et alJl2006h. obser- 
vations from the ultraviolet fluxes (IHicks & Mushqtzkvl2005t : 
Ifficks et alj [ 2009t iPipino et al.ll2009bt IWang et alj|2009t) . the 
infrared fluxes ( Egami et al.ll2006t lO'Dea et alj 120081). opti- 
cal photometry Bildfell et al.ll2008t iRaffertv et al.ll2008h and 
optica l line emission (ICrawford et all 119991: lEdwards et alj 
120071) of the BCGs suggest that recent star formation is still 
taking place in some of these objects. In particular, about 
one third of the X-ray selected BCGs associated with cool- 
core clusters exhibit recent star formation. This indicates that 
AGN feedback may not completely suppress star formation 
in BCGs at the center of cool-core clusters. Here we assume 
that 10-20% of stars may form from cooling gas at cluster cen- 
ters. 7 Using a — 4 and / con =0.1 for the present-day BCGs, 
we expect that the UV flux would be enhanced by 40% due 
to the He sedimentation effect (from eq. [HI). The FUV - V 
color, therefore, is ~ 0.4 mag bluer. 

Therefore, one of the predictions of our model is that the ef- 
fect of He sedimentation, and hence the U V flux, is the largest 
for the BCGs in cool-core clusters than in non-BCG elliptical 
galaxies. Systematic studies of the UV flux strength in BCG 
and non-BCG elliptical galaxies can test our models. 

On smaller scales, recent observations suggest that some 
globular clusters (GCs) in our Galaxy, such as a) Cen, re- 
quires He-enhanced population (Y = 0.4) to explain both the 
blue main-se quence population as well as the ho t horizontal 
branch stars dLee et al.l l2005b; Pio tto et alj 12005b . However, 



we point out that the global He sedimentation process is likely 
unimportant for GCs. It is because that GCs pass through the 
Galactic disk on a timescale of ~ 0.1 Gyr and each cluster- 
disk interaction would wash out the accumulated gas in the 
GCs via ram pressure stripping. Since the interaction interval 
is much less than the He sedimentation timescale (~ 1 Gyr for 
a 10 6 M o GC), there is not sufficient time for He sedimentation 
to build up significantly high He abundance. 

4. CONCLUSIONS AND DISCUSSIONS 

In this work, we propose that the He sedimentation can 
significantly enhance the UV signal observed in BCGs. We 
show that the He sedimentation process can produce signifi- 
cant amount of He abundance (AY > 0.25) in central regions 
of hot, massive galaxy clusters. Stellar populations in BCGs 
may contain UV bright EHB stars when they are old. The 
model predicts that the UV upturn phenomena increases with 
X-ray temperature or mass of galaxy clusters. 

Our sedimentation model makes several unique predica- 
tions that are testable with current GALEX+SDSS observa- 
tions. 

1 . The model predicts a strong correlation between the UV 
upturn strength of BCGs and the X-ray temperature of 
host galaxy clusters. 

2. The model also predicts that the UV upturn phenomena 
should be more pronounced in BCGs, compared to non- 
BCG elliptical galaxies. 

3. Some correlations are expected between the UV upturn 
strength of BCGs and the dynamical state of the BCGs 
and their host galaxy clusters. Major mergers, for exam- 
ple, can destroy He-rich cluster cores, producing some 
scatter in the UV upturn among BCGs. In this case, we 
expect a correlation between the size of UV upturn and 
dynamical state of clusters (e.g., measured using X-ray 
morphology or multiple X-ray peaks). This is another 
expected observational signature of our model. 

Analyses of a large sample of galaxies from the 
GALEX+SDSS along with average temperatures of their host 
X-ray clusters can provide important test of this model and 
may shed new light on the origin of the UV upturn phenom- 
ena as well as the physics of the formation and evolution of 
massive galaxies in the universe. 
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